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Abstract. For some blazars, the gamma-ray absorption features due to pair-production on the 
Extragalactic Background Light (EBL) are fainter than expected. The present work reviews 
the main models that could explain this paradox, with emphasis on conservative ones, that 
do not include any new physics. The models that are intrinsic to the source, do allow a 
very hard primary spectrum, but fail to explain a regular redshift dependence of the anomaly 
starting energy. The model that includes a contribution from secondary photons produced by 
cosmic rays (CR) near the Earth seems to require a well collimated CR beam, what is hard to 
achieve. Finally, the model with secondary photons produced in electromagnetic (EM) cascades 
initiated by primary gamma-rays is considered. In principle, it allows to decrease the statistical 
significance of the anomaly and, while requiring quite low EGMF strength B, does not contradict 
to most contemporary constraints on the B value. Additionally, it is shown that the recently 
observed correlation between directions to hard gamma-ray sources and voids in the Large Scale 
Structure is a natural feature of the EM cascade model. 


1. Introduction 

The number of detected extragalactic gamma-ray sources has greatly increased during the 
last decade. Blazars — gamma-ray loud active galactic nuclei (AGN)— constitute the great 
majority of these sources. By the start of 2015, ~50 blazars were discovered by ground-based 
detectors [T] (these instruments typically work at very high energy (VHE) range E >100 GeV), 
and ~ 10^ blazars were observed by the Fermi LAT gamma-ray space telescope at E= 100 
MeF-100 GeV [2]. Besides investigation of AGN intrinsic properties, these observations allow 
to study extragalactic gamma-ray propagation as well. High energy photons are subject to 
the 77 —>■ e'’'e“ process on the Extragalactic Background Light (EBL) photons [3]-[l] with 
subsequent formation of electromagnetic (EM) cascades. Thus, primary gamma-rays with 
E > Eabs {Eabs = E{T-y-y = 1); Xy-j, is the optical depth of the pair-production process) are 
effectively “absorbed”, and their energy is transferred to electrons, positrons (in what follows 
called simply “electrons”) and photons with comparatively low energy. Absorption features of 
the 77 —>■ e'’'e“ process were already observed by the Fermi LAT instrument [5], as well as the 
H. E. S. S. Gherenkov telescope [6] with high statistical significance: ~6(T and 8.8(T, respectively. 

However, it appears that for some blazars modification of the gamma-ray spectrum in the 
optically thick regime has a somewhat anomalous character [7]. By extrapolating the spectrum 


from the optically thin regime, [7] found that the distribution of the flux points scatter around 
the predicted intensity is different for the 1< <2 and >2 regions, and a strong indication 

(significance Za = 4.2ct) that the observed absorption is smaller than the assumed one was 
obtained. This effect, observed at VHE range, was called “the pair-production anomaly” by [7]. 

In fact, an anomaly of such kind is a long-standing problem of VHE gamma-astronomy [ 8 ]- 
[9]; a number of solutions was proposed, including quite exotic ones, such as violation of Lorentz 
invariance (LIV) [TU]-[TT].p]. and oscillations of photons to Axion Like Particles (ALPs) (e.g., 
m-M)- A search for new physical phenomena of such kind would require that all possible 
background processes are well known and appropriately accounted for. One such process is 
emission of secondary photons in EM cascades. On a subsample of spectra analysed in m, it 
was recently shown that inclusion of this secondary component to the spectral fitting procedure 
indeed allows to considerably decrease the statistical significance of the anomaly [K]. In the 
present work the cascade model of the VHE anomaly in blazar spectra is compared with other 
models of the anomaly. The paper is organized as follows: section 2 contains a brief overview of 
possible processes that might relax the anomaly; the EM cascade (EMC) model is described in 
section 3; section 4 reviews some contemporary constraints on the Extragalactic Magnetic Field 
(EGMF) strength. Finally, some conclusions are presented in section 5. 

2. The models of the VHE anomaly 

2.1. EBL intensity models 

The simplest assumption that could relax the anomaly would be that the EBL energy density is 
overestimated. If this is the case, the intrinsic spectrum would experience weaker modification 
than is believed, and for sufficiently low EBL density the VHE anomaly would be relaxed (note, 
however, that the dependence of the anomaly significance on the EBL intensity may have a 
non-trivial behaviour |16|). 

We are mostly interested in the wavelength region A=l“10 fim that roughly corresponds to 
1~10 TeV minimal energy of high-energy photon that is subject to the 77 —^ e"*“e“ process. 
Many EBL models exist, including m-m-. in the range of A=l“10 //m the model of that 
was utilized in [7j, predicts one of the lowest EBL intensity among [I7]-[53]. The [T7] model is, 
in fact, a sort of lower limit on the EBL intensity. A well justified conclusion that such a limit 
must be revised would likely require a large amount of dedicated observational work, therefore 
we move to other possible explanations of the anomaly. 

2.2. A list of effects that could relax the anomaly 

Besides the EBL properties, a number of factors exist that may have some relation to the nature 
of the VHE anomaly. A non-exhaustive list of them is given below. 

1. Effects, connected to the source. 

1.1. Intrinsic spectrum may have a pile-up resulting from inverse Compton (IC) scattering on 
an ultrarelativistic outflow, thus increasing observed intensity at high energies [25]. 

1.2. Internal absorption of primary photons on radiation field near the source, under certain 
assumptions, again, may produce a very hard intrinsic spectrum [26| with a pile-up. 

1.3. Another known mechanism of a hard spectrum formation is synchrotron radiation of 
electrons produced by ultrahigh energy (UHE) protons near the source [77] . 

1.4. It is believed that blazars have well collimated jets pointed towards the observer (the scheme 
of a typical blazar geometry taken from |55| is drawn in figure 1). Magnetic field in a typical jet 
likely has some turbulent component, that may scatter charged energetic particles (electrons and 
hadrons). This scattering angle is usually larger for charged particles with comparatively low 
energy, that emit predominantly low-energy photons. Thus, the beaming factor for high-energy 
photons would be larger, again increasing observed intensity at high energies. To the author’s 
knowledge, this effect was never discussed in connection to the VHE anomaly yet. 
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Figure 1. A sketch of a typical AGN geometry [28]. 1 — 
central compact object (black hole with accretion disk), 2 — 
dusty torus around the central object, 3 — well collimated 
relativistic jet with regular magnetic field Bj^t that can have 
some turbulent component 5Bjet, 4 — a “lobe” at the end of 
the jet that likely contains a highly turbulent magnetic field 
region. 


However, it was observed that the energy at which the anomaly shows up, exhibits strong, 
regular dependence on redshift (this energy is defined by the relation Ea ~ E{t^^ = 2) [7]). The 
effects, intrinsic to the source, are not expected to display such a dependence, thus they can not 
constitute the only physical reason of the anomaly. 

II. Propagation effects. 

11.1. LIV effects, in principle, may suppress the 77 —>■ e^e“ process [ID] . This mechanism is 
expected to operate above a certain fixed energy, and so it is disfavoured for the same reason as 
the source-induced effects (see discussion in my 

11. 2 . 7 -ALP oscillations (and vice versa) in magnetic field is another mechanism that makes 
the primary gamma-rays able to avoid strong absorption m- Even if we assume the existence 
of ALP and that the 7 -ALP oscillation process is indeed allowed, there is another potential 
difficulty: the efficiency of direct (7 —)• ALP) and reverse (ALP —> 7 ) conversion is required 
to be high enough for many sources, independently of direction to a source, and, therefore, of 
the Galactic magnetic held parameters in this direction, as well as of the source’ magnetic held. 
The recent data didn’t allow the authors of |29] to put any constraints on the ALP mechanism, 
but the same work shows that such constraints are obtainable in the future. 

11.3. If blazars accelerate hadrons as well as electrons, secondary gamma-rays produced by the 
former may contribute to observed spectrum and thus relax the anomaly m-m- This model 
is discussed below. 

11.4 As was already mentioned, the most basic process that may inhuence the shape of observed 
spectrum is emission of secondary photons in electromagnetic cascades | 8 ], [25], [32], [15]. This 
process is discussed in section 3. 

2.3. The CR beam model 

Secondary photons produced by cosmic ray beam (by means of pair-production and 
photohadronic processes with subsequent development of electromagnetic cascades) relatively 
near to the Earth may enhance the intensity in the optically thick regime. Secondary gamma- 
ray spectra, typical for this mechanism, were studied in [33]-|36|. The case of a very distant 
source (with redshift z ~1) was considered in [37]. Time structure of VHE gamma-ray signal was 
discussed in [33] • GR beam may also produce an observable amount of neutrinos [33] -[34].[39]. 

This model of the VHE anomaly has many attractive features: it allows to explain the 
observed spectrum without invoking any exotic processes such as 7 -ALP oscillations. Moreover, 



the shape of the secondary photon component, assuming sufficiently low strength of the EGMF 
B < 10“^^ G, depends mainly on the overall normalization of the primary CR spectrum. 
However, the model is not out of difficulties. As an example, let us consider the lES 1101-232 
blazar spectrum (redshift z=0.186). Figure 2 shows a spectral energy distribution (SED) of the 
source measured by the H. E. S. S. Cherenkov telescope (red circles) together with uncertainties 
(dashed red lines); the assumed intrinsic spectrum is shown by black curve, the fit to the 
observed spectrum is shown by green line (for detals of calculations, that were performed with 
the publicly-available code ELMAG 2.02 [40j and the EBL model of |T7], see |15jl. The intrinsic 
spectrum curve here is normalized to the absorbed one at E= 200 GeV, where absorption effects 
are already small. The fit here doesn’t include any secondary photon component. 



Figure 2. The fit to the lES 1101-232 
SED [H] without account of any secondary 
emission [15] . 



Figure 3. The fit to the lES 1101-232 SED 
including cascade component (blue curve). 
Figure from |15] . 


To explain the data in the framework of the CR beam model, an additional component from 
secondary photons is needed (the corresponding region of the spectrum is denoted by thick blue 
dashed lines). Now let us put a simple lower limit to the ratio of total power of accelerated 
CR WcR to the same quantity for VHE electrons WvHE-e = WvHE-'y/fiCi where fjc is the 
fraction of energy of VHE electrons transferred to VHE photons by means of IC scattering. An 
important note is that gamma-rays produced by electrons deep inside blazar jets are expected 
to be highly beamed. On the other hand, a region of highly turbulent magnetic field (a “lobe”) 
is expected to exist at the end of a jet, and CR particles are likely to be isotropized before they 
could produce secondary photons close to the Earth. Indeed, such “lobes” are observed in a 
number of radiogalaxies, and are expected to be able to scatter (or even confine) protons up to 
10^^ eV [l2]. Therefore, if the observer is located near to the jet’s axis, the observed flux from 
CR is effectively deamplified to a factor of fbeam ~ assuming the angular radius of the 

jet 0jet=l-3°. Additionally, only a small fraction fint ~ 10“^ [3T] of CR energy is converted to 
photons near the observer. 

Assuming the intrinsic VHE gamma-ray spectrum dN/dE ~ E~‘^ (what is close to the shape 
of black curve in Figure 2) from 100 GeV to 100 TeV, one may estimate that the secondary 
photons constitute a fraction fg^c ~ 10“^ of energy of intrinsic VHE photons. Therefore, 
assuming fjc = 0.1, we find that Wcr/Wvhe-s ~ fsec ■ fbeam ■ ific/fint)= lO^-lO"^, depending 
of the fbeam value, i.e. acceleration of hadrons is very effective, while for VHE electrons it 
is not. This challenges conventional models of gamma-ray blazars. Additionally, assuming 







WvHE-e = 10 "^^- 10 ^®er( 7 /s, we find Wcr= W^^erg/s, a very large value. The last values 

are similar to those obtained by [43| . 

3. The EMC model 

3.1. General considerations 

The primary photon absorption process 77 —>■ e'*"e“ has threshold energy E.yo = rUg/e, where rUg 
is electron mass, and e is EBL photon energy (we assume c=l here). Most of pair-production acts 
occur not far from the threshold; in this case each produced electron receives energy Eg ~ E.yQ/2. 
Electrons, in their turn, produce secondary (cascade) photons by IC scattering. This process, 
in contrast to pair-production, doesn’t have kinematical threshold, and occurs mainly on dense, 
low-energy cosmic microwave background (CMB) photon field; the energy of secondary photons 
E-y-casc ~ 3 • {Eejl TeV)"^ GeV [iQ]. In this work we consider the case 100 GeV < E.yo < 100 
TeV, thus E.y-casc < 7-5 TeV « E^q-^ for the case of E^q = 10 TeV E^^^casc = 75 GeV, again 
<< E.yQ. Therefore, the cascade process occurs in a peculiar, “degenerate” regime, and a typical 
number of cascade generations is small, njnt = 1 - 2 . 

In absence of EGME the angular radius of cascade Oint ~ ^/nt/7") where T = Eelroe > 10^ 

— Lorentz factor of a cascade electron with Ef, >100 GeV = 10® MeV. As we have seen, 
9Jet ~ 1 ° ~ 10 “^ rad; then Ojnt < 10 “® << Ojet, therefore the cascade process weakly modifies 
the angular distribution of photons from the source and ID assumption that we shall use is well 
justified. Recently it was argued that plasma beam instabilities could dominate the energy losses 
for intergalactic cascade electrons |H]. However, this mechanism is under debate; it was noted 
that the relativistic length contraction effect “compresses the electric and associated magnetic 
fields” around charged particles, and so, in fact, they do not constitute a plasma [35]. For this 
reason, we do not include any electron energy loss process, except IC, in our calculations. 

3.2. The signatures of cascade emission 

Now let us illustrate the basic features of the EMC model of the VHE anomaly (figure 3). A fit to 
the measured SED was performed m. including the secondary component from electromagnetic 
cascades. The notations are the same as in figure 2, but the total model SED is now denoted 
by blue curve; the model SED without account of the cascade component is drawn by green 
curve. The impact of cascades to the SED is the difference between the blue and green curves. 
Blue vertical line denotes the “gamma-ray horizon” (the energy, where Tyy=l)', red vertical line 

— the energy, defined by condition T.y.y=2. As in figure 2, the intrinsic SED in figure 3 was 
normalized to the absorbed one at E= 200 GeV. 

The above-discussed rule E^^casc « is clearly visible in figure 3, and the cascade 
component contributes mainly in the optically thin regime. However, to achieve a good fit, we 
were compelled to take a very hard primary spectrum, much harder, than in figure 2 for the 
same object neglecting the secondary component. Therefore, the model intensity in the optically 
thick regime in figure 3 (with account of the cascade component) is much higher, than in figure 2 
(without this component), and for >2 the fit is in much better agreement with observations. 
It is quite remarkable that an exceptionally good fit of the shape of the spectrum in figure 
3 (13 energy bins) was achieved with a very simple parametrization of the primary spectrum 
dN/dE (X Eq^ ■exp(—Eo/Ec} with only two parameters 7 and Ec. As well, the significance of the 
anomaly for >2 drops from 2.1 cr for the fit shown in figure 2 (without secondary component) 
to only 0.5 cr for the new fit that includes the cascade component (for the significance evaluation 
method see [T^l. 

A possible influence of ECMF to the observed spectrum is schematically shown in Figure 3 
by red dashed curve. In order to not spoil the fit, the EGME strength B should be sufficiently 
small, below some value Bmux- Recently A. Neronov et al. |l 6 ] performed calculations of this 
“magnetic cutoff” shape for the case of blazar Mkn 501 {z= 0.034), and found that for B = 10“^® 


G the cutoff starts to be detectable below E ~ 100 GeV. Given that, as was noted above, the 
number of generations in cascade is small, and that the last low-energy point of the observed 
spectrum, located sX, E ^ 200 GeV, doesn’t show a prominent magnetic cutoff, we can estimate 
Bmgx ~ (1 — 2) • 10“^® G. Plausibility of such B values will be discussed below in section 4. 

A spectrum of blazar in the considered version of the EMC model has three characteristic 
features, or signatures: a) low-energy magnetic cutoff b) the “dip”, or a kind of an “ankle” 
that is usually located in the energy region where 1 < < 2 c) high-energy cutoff (not shown 

in figure 3). The third feature, the high-energy cutoff, was already studied in [33], [IZ]- The 
detailed properties of this signature are not known, but the general conclusion that, in principle, 
allows to discriminate between the CR beam model and the EMC model, is the following: for 
large (say, >2-3) the high-energy cutoff is much more marked for the case of EMC model 
ESI, m- The dip or ankle feature, as well, exists in both the CR beam (see, e.g., [33]) and 
the EMC model. However, the low-energy magnetic cutoff is much more marked in the EMC 
model when compared to the CR beam model, unless the secondary component produced by CR 
dominates the entire VHE spectrum. We will see in section 4 that existing hints at intergalactic 
cascades tend to favour the EMC model. 

4. Some constraints on the EGMF strength 

As we have seen, by now there exist only two conservarive models that could relax the VHE 
anomaly in blazar spectra, not to mention the option that the EBL intensity models must be 
revised: the CR beam model and the EMC model. Both require sufficiently low EGME strength 
in voids to be plausible: B < 10“^^ G for the former [38|, and even lower values (usually around 
B < 10“^® G) for the latter. 

Any conclusive measurement of EGME strength is still absent; some constraints on this 
quantity are presented in figure 4 (usually they are obtained for the EGME coherence length 1 
Mpc [3S|)- Upper bounds on the B value are scarce and rather weak; [33| obtained a constraint 
B < 10“® G. B = 2 ■ 10“^^ G was found to be sufficient to explain magnetic fields in galaxy 
clusters [3D|. As well, we present a number of results obtained from non-observation of cascade 
component in blazar spectra at comparatively low energies, E <100 GeV. These constraints 
are highly model-dependent; m even found that the zero EGMF hypothesis cannot be firmly 
rejected. Anyway, we present the cascade constraints obtained by [52|-[53] in figure 4; these 
lower limits are in the range B > 10“^®-10“^^ G. By analysing the angular distribution of 
arriving photons, m performed a search of a magnetically broadened cascade pattern and was 
able to exclude the range of H = 3 • 10“^®-10“^^ G at 99 % C.L. Finally, by thick red horizontal 
line we show the B = 10“^^ G value, above which both the CR beam model and the EMC 
model are ruled out [38|. Future instruments with high sensitivity (e.g. CTA [56]), or with good 
angular resolution, such as emulsion gamma-ray telescope m or the GAMMA-400 detector |58| , 
would be able to employ the same method as [55| to put more tight constraints on the EGMF 
strength. Other EGMF searches based on stacking analysis of sources’ angular distribution [59] 
or diffuse gamma-ray sky studies |60j- m do exist. 

Finally, let us mention some recent results that actually give hints at the EM cascade presence 
in blazar spectra. Analysing the Mkn 501 spectrum, [36] found that it could be explained by 
the intergalactic cascade hypothesis with B = 10“^^-10“^® G. Very recently, it was shown 
that gamma-ray sources with hard spectra are predominantly located in directions to voids [62] . 
Moreover, it was found by [62| that the EBL intensity fluctuations are not sufficient to explain 
the observed effect. Therefore, as first suggested in m, the most likely explanation is that 
|62j actually observed the low-energy magnetic cutoff of secondary component, produced by EM 
cascade. This feature of the EMC model is schematically shown by dashed red line in figure 
3. For the case of other objects, that do not point to voids, the mean B value is likely much 
greater, and the cascade component may be greatly suppressed. For this reason, these latter 



Figure 4. Some constraints on the EGMF 
strength in voids, n denotes the number of 
corresponding work. n=l — |49] : n=2 — |50j : 
n=3 — [52]; n=4 — [53]; n=5 — [5l], the case 
of the EBL model m-, n=6 — [5l], the case of 
the EBL lower limit from direct source counts; 
n=7 — [55] n=8 — [38]. 


objects do not show the spectral hardening observed by [U2] . 

In principle, two other explanations of this effect are possible. The first is that the 
electromagnetic cascade upscatters the secondary photons to high enough energy to produce 
observable hardening. However, this is highly unlikely, given that a typical energy of secondary 
photon is much smaller than of the primary one, as discussed in subsection 3.1. The other 
explanation, that the effect observed by [62] is caused by the CR-initiated cascade, is unlikely 
as well, unless the secondary component from these cascades dominates the VHE spectrum. 

5. Conclusions 

In the present work a brief review of conservative models that could explain the apparent 
faintness of the pair-production features in blazar spectra, was given. No existing model was 
found to be completely free of difficulties. The source-intrinsic models, while they are able to 
form a very hard intrinsic spectrum, do not explain other features of the anomaly (subsection 

2.2) . The model with production of secondaries by CR near to the Earth seems to require either 
a very high Wcr/Wyhe-c ratio, or extremely good collimation of accelerated nuclei (subsection 

2.3) . Finally, the EM cascade (EMC) model was considered (section 3), that naturally explains 
some spectra with only two free parameters, doesn’t contradict to contemporary constraints on 
the EGMF strength, and predicts the magnetic cutoff feature that may be already observed [62]. 
Future observations will help to test this model. 
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